EXAMPLE ON THE APPLICATION OF FORCE AND 
DISPLACEMENT METHODS ON THE ANALYSIS OF STATICALLY 
INDETERMINATE FRAMES 


1.0 Introduction to methods of structural analysis 


It is a well known fact that there are two basic methods of structural analysis of indeterminate 
structures which are; 


1. Flexibility methods (also known as force methods, compatibility methods or the methods of 
consistent deformation) and 
2. Displacement methods (also known as stiffness or equilibrium method) 


Each method involves the combination of a particular solution which is obtained by making the 


structure statically determinate, and a complementary Sotuft?/ y which the effects of each individual 
modification is assessed. In the force methods, the behaviour of the structure is considered in terms of 
unknown forces, while in the stiffness method, the behaviour of the structure 1s considered in terms of 
unknown displacements. By implication, both analysis methods always involve reducing the structure 
to a basic system (a determinate system). In the force method, the basic system involves the removal 
of redundant forces, while the stiffness method involves restraining the joints of the structure against 
displacement. 


Whenever we are using the force method, our basic system is dependent on the degree of static 
indeterminacy, which a system that must be statically determinate and stable. The choice of redundant 
constraints to remove is based on experience, geometrical configuration, stability of the structure, and 
ease of analysis. Degree of static indeterminacy can be computed using; 


D = (3M + R) - 3N --------------------------- (1) 
Where; 
D = Degree of static indeterminacy 
M - Number of members 
R = Number of constraints (reactions) 
N = Number of nodes 
Alternatively, the following relation can be used; 

D = R = 3 — S -------------------------- (2) 

Where; 
D = Degree of static indeterminacy 
R = Number of constraints (reactions) 
S = Number of special conditions (such as an internal hinge) 


After selecting a good basic system, we now replace the redundant supports with unit loads and 
evaluate them one after another. In this step, we calculate the deflection corresponding to each 
redundant force separately due to applied loading and other redundant forces from force-displacement 
relations. Deflection due to redundant force cannot be obtained without knowing the magnitude of the 
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redundant force. Therefore, we apply a unit load in the direction of redundant force and determine the 
corresponding deflection. Since the principle of superposition is valid in elastic analysis, the 
deflections due to redundant force can be obtained by multiplying the unknown redundant with the 
deflection obtained from applying unit value of force. Now, we calculate the total deflection due to 
the applied loading and the redundant force by applying the principle of superposition which must be 
compatible with the existing boundary condition. 


For more than one set of redundant forces, we construct a set of simultaneous equations with the 
redundant forces as unknowns and flexibility coefficients as coefficients of the equations. These 
flexibility coefficients are also called the influence coefficients. The total number of equations equals 
the number of unknown redundant forces. These set of simultaneous equations are called the 
canonical equations, and they are of the form shown below; 


O11 X4 EE 015X5 Taat Ojo X, + AP =0 
051 X4 RE 055 X5 potias OonXn + AP =0 


Aa X4 T Òn2 X» WoT Dos a A,P =0 
Where; 
0;- deformation at point 7 due to a unit load at point j 
Xi = Redundant force at point i 
A;P = Deformation at point į due to externally applied load 


Mohr's integral 5, - f, ns (dS ecann (4) 


Where; 

0, - deformation at point 7 due to a unit load at point j 

M - Moment due to externally applied load or another redundant load 
m = Moment due to unit load (virtual load) 


Displacement method 


Analysis of statically indeterminate structures by the displacement method in the canonical form 
begins with determining the degree of kinematical indeterminacy. The fundamental approach in the 
displacement method is the opposite of force method. First of all, we calculate the deformations at the 
ends of the members and then the internal forces in the members. Thus, the primary unknowns in the 
displacement method are the displacements. 


Analysis of any statically indeterminate structure by the displacement method begins with 
determining the degree of kinematical indeterminacy. Generally, the degree of kinematical 
indeterminacy n of a structure is determined by the formula; 


Where; 


n, is the number of unknown angles of rotation of the rigid joints of a structure, and 
nq is the number of independent linear displacements of the joints. 


In general, the degrees of kinematical and static indeterminacy are not equal. The primary system of 
the displacement method is obtained from the given one by introducing additional constraints to 
prevent rotation of all rigid joints and all independent displacements of various joints. These 
introduced constraints are shown by the shaded squares and the double lines, respectively. 
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Primary unknowns Zi (i = 1, 2, 3, ... n) represent displacements of introduced constraints (angles of rotations 
and/or linear displacements of various joints of a frame. The number of primary unknowns, n, for each structure 
equals to the degrees of its kinematical indeterminacy. 


The canonical equations of the displacement method will be written as follows: 


kiiZi or ki2Z2 Wase T kipZn + KP =0 
koiZ4 + ko5Z5 "ust kona T KP =0 


kwZ T kwLo PT KanZn T K,P =0 
Where; 
kj; - Unit reaction (force or moment) at point 7 due to a unit displacement at point j 
Xi = Real unknown displacement at point 7 
A;P = Reaction at point i due to externally applied load 


3.0 Guide to choice of method 


It is expected that the reader should by now have known the basic principles underlying both analysis 
methods. Often times in the classroom, the lecturer often specify the method he wants the students to 
use, based on the scheme of work. Apart from that, sometimes it is required that manual methods be 
used to verify some results produced by computer programs, or when you are verifying the results 
produced using a MATLAB code you have written. The following are simple guides to choice of 
method. 


1. Check the degree of static and kinematic indeterminacy using equations | and 5 above. The 
one that has the less value is usually the most favourable when you are using your classroom 
non-programmable calculator. However, when you have programs that can help you solve the 
matrices (e.g. Microsoft excel or MATLAB) such advantage could be written off. 

2. Check the loading configuration on the structure. The more complex the loading, the less 
advantageous force method becomes over the displacement method. This is because in the 
force method, the graphical method (use of moment diagrams) or the integration method 
(Mohr’s integral) is usually used. The more complex the loading, the more complex the 
diagrams are expected to be, and of course the more number of sections you are expected to 
cut in order to generate moment equations for your integration. But if the displacement 
method is used, all loads are converted to FIXED END MOMENTS and SUPPORT 
REACTIONS, which is quite handier than the former. 

3. When you are seeking to obtain the rotation at the nodes, the displacement method gains more 
advantage over the force method because the solution of the canonical equation gives the 
displacement of the structure. However, if force method is used, it is imperative to obtain the 
final bending moment diagram of the structure, before proceeding. This is also the case when 
you are looking to obtain the side sway of the frame. 

4. Based on the nature of the analysis methods, it is faster to obtain support reactions of frames 
using force method since the solution of the canonical equations give the solution. However 
in displacement method, support reactions can be obtained from 1* principle after analysis of 
the bending moments. A variation of this is in the direct stiffness method, which is better 
handled using a standard computer program. It is not handy for manual calculations at all. 

5. The force method handles frames of non-symmetric configuration better than the 
displacement method when computing manually. For instance, it is easier to handle frames 
that are made up of slanted members with force method than with displacement method since 





Page 3 


nothing changes about the way moment diagrams are computed. When slanted or non- 
symmetric members are involved in displacement method, the side sway mechanism can get 
complicated, and errors can creep in. 

6. When writing structural analysis programs in computer, displacement method (specifically 
direct stiffness method) is beyond the competition. This simply requires finding a way of 
obtaining and handling the structure's general stiffness matrix because while if force method 
is employed, there is a big task of programming how to reduce the structure to basic system 
(which is a system that is statically determinate and stable). For some structural 
configurations, it can be a quite daunting task. Even during manual analysis, it is common to 
see some students getting the basic system wrong. But the basic system in displacement 
method is easy to conceive. 


4.0 Analysis Examples 


Problem 1.0: A frame is loaded and supported as shown below, obtain and draw the internal stresses 
diagram (moment, shear and axial) for the externally applied load. Also obtain the side sway of the 
frame using; 


a) Force method 
b) Displacement method 


(EI = Constant) 


10 kN 








| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
X. 200m y 200m J 4.00m X 
| 


Solution 
NB: DIAGRAMS ARE NOT TO SCALE. READ ASSIGNED VALUES 
(a) By force method; 
Step 1: Determine the degree of static indeterminacy. 


D-R-E-s where; 
D = Degree of indeterminacy 
R = number of support reactions in the system = 5 
E = number of equations of equilibrium = 3 
S = any special condition, e.g. internal hinge = 0 
Therefore, N25-3-0z2 
Therefore, the structure is indeterminate to the 2" degree 





Page 4 


And we are expected to remove two redundant supports, and solve a simultaneous equation containing 
two unknowns in our canonical equation. 


Step 2: Reduce the structure to a basic system 
As we realised from the first step, we will need to remove two redundant supports in order to make 


the structure statically determinate. However, we must ensure that the selected system must be stable. 
If we remove the roller at supports E and F, we will obtain a good basic system. So in this example, 
we are removing the reactive forces at supports E and F and replacing them X; and X» which will be 


assigned unit values. The basic system is shown below. 





Step 3: Analysis of the various load cases (X,, X; and externally applied load) 


The values X, and X, represent the redundant forces that have been removed from the system, and 
they will be assigned unit values in order for us to progress in our analysis. We are going to treat each 
of them as independent load cases on the structure. Our main interest in this discuss will be on the 
moment diagrams that they produce. You are expected to be familiar with the analysis of statically 


determinate frames by now. 


Case 1 (X, = 1.0, X, = 0) 


6.00m 





y 4.00m N 4.00m y 
| | 
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Case 2 (X, = 0, X, = 1.0) 














" 


i 
, 


8.0 


6.00m 


UAI 






8.0 





Loading the basic system with externally applied load 


20 kN 






9 KN/m 
Ey vyyyvyvyvyvevvyvvvv 





5 2.00m y 2.00m N 4.00m EN 


The resulting bending moment diagram due to external load is given below; 


160 kNm 


3.00m 


3.00m 





190 kNm 2.00m 2.00m 4.00m 
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Step 4: Computation of influence coefficients 


044 = Deflection at point 1 due to unit load at point 1 


4.0 4.0 





à - (ix 4x 4x 4) + (4x 4 x 6) = 117.333 


612 =012 Deflection at point 1 due to unit load at point 2 





512 = (4x 8 x 6) += x [AX (2x 8) + 4) x 4] 
547 = 84 = 192 + 53.333 = 245.333 


622 = Deflection at point 2 due to unit load at point 2 





85, = (=x 8 x 8x 8) + (8x 8 x 6) = 554.667 


AP: Deformation at point 1 due to externally applied load 


2.0 


2.0 





AiP = —= x 2 [160((2 x 4) + 2) + 80((2 x 2) + 4)] - 2x [2 x ((2 x 80) + 40) x 2] — (40x 


160 x 3) —= x [4 x (160 + 190) x 3] = — 4900 
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^P: Deformation at point 2 due to externally applied load 
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A;P - - (2x 40 x 4x 4) - 2x2 [160(2 x 8) 6) + 80(2 x 6) + 8)] -+ x 2 [80(Q2 x 


6) * 4) -- 40((2 x 4) + 6)| — (8 x 160 x 3) —= x [8 x (160 + 190) x 3] =— 10520 


Step 5: Canonical equations 
The appropriate canonical equation is 


014X, + 042X204 AP =O 
091X1 RS 025 X» + A P = 0 


117.333X, + 245.333X, = 4900 
245 .333X, + 554.667X, = 10520 
Solving simultaneously, 
X, = 27.995 KN 
X» = 6.584 KN 


Note now that the values X, and X; represent the actual vertical support reactions occurring at roller 
supports E and F. Since we have known these values, the structure has become determinate, and we 
can analyse it by the law of statics, or progress using the force method. In this example, we are using 
the force method to obtain the final moment diagram, and not the shear and axial forces (see further 


examples for shear and axial). 
Step 6: Final internal stresses (moment only in this example) 


M; = M,X, + MX; + Mo 

Ma = (27.995 x 4) + (6.584 x 8) — 190 = —25.348 KN.m 
Mg- Mc? = (27.995 x 4) + (6.584 x 8) — 160 = 4.652 KN.m 
Mc“ = (27.995 x 4) + (6.584 x 8) — 160 = 4.652 KN.m 

Mp = (27.995 x 2) + (6.584 X 6) — 80 = 15.494 KN.m 

M; = (27.995 x 0) + (6.584 x 4) — 40 = —13.664 KN.m 
Mr= 0 


Determine the maximum moment on span E-F using any method of your choice 





Step 7: Plot the final internal stresses diagram (moment only in this example) 


13.664 kNm 





15.494 kNm 






4.652 kNm 





To analyse the same frame by displacement method, we make the structure kinematically determinate. 
To start with, we consider the number of unknown nodal displacements that the structure can undergo. 
A little consideration will show that the structure is kinematically indeterminate to the 4" degree. 
These unknown displacements are the rotations at nodes C, E, G, and the lateral translations at F (side 
sway). To make the structure kinematically determinate, we rigidly fix up nodes C (labelled 1), E 
(labelled 2), and F (labelled 3), and also put a constraint at nodes F (labelled 4) to prevent lateral 


displacement. This is shown in the figure below; 


Z4 Zə Z3 


e —s 





6.00m 
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Analysis of case 1; Z; = 1.0 


Applying a unit rotation at point 1 


— Ka 





y 4.00m ` 4.00m S 
Ky, = += = L667EL Ka == =0.5E; Ka =0; Ky == = — 0.167EL 


Analysis of case 2; Zz = 1.0 


Applying a unit rotation at point 2 


—— Ka 





y 4.00m ` 4.00m ` 
K= = = 0.5EL Kə» = T + = =1.75EL Ky =0:Ky =0 


Analysis of case 3; Z; = 1.0 


Applying a unit rotation at point 3 
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Ki = 0; K33 = — = 0.5EI; K33 = T = EL K;5 = () 


Analysis of case 4; Z4 = 1.0 


Applying a unit lateral displacement at point 4 


6.00m 





Ki => = — 0.167EE; Ka, = 0; Ku, = 0; Ku, = = = 0.0555EI 


External load 





P PL 
kpo = >= 
-.PL ql’ 20x4 5x4 
Kop= 8 8 8 79 
K35o— 0 
P 10 
Kz (reaction from horizontal load) = — iio —5 KN 
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The general canonical equation is given by; 


kZ + KioZo +...4+ KinZn + Kip = O 
ko1Z4 + ko5Z5 east ko»Zn + Kop = 0 





; E E NES. Do --------------------- (10) 
KniZi T Kn2Zo "RT Kinla + Kap =0 
Hence, 
1.667 05 0 -—0. ~ Ly 2.9 
0.5 175 0 Z2|_| 0 
0 0.5 1 Z3 0 
—0.167 0 0 0. an ZA 5 
On solving; 
Z= Lo radians; Z» = A: radians; Z4 — E radians; Z4 — —— metres; 
5.7 Final Moment Values: 
Mana = Mo + M,Z) + M>Z> + M3Z3+ M,Z, aaa a aa (10) 


17.0930 /2EI 141.2791 /6EI B 
Ma = —7.5 4 ZE (= )+0+0 — 1412791 (SS) = — 25.348 KNm 














B. 17.0930 4EI B 141.2791 /6EI 2. 
Mo” = +7.5 += (—=)+0+0 = (>) = 4.651 KNm 
R 17.0930 /4EI 4.8837 /2ElI 
Mc = —10 + 222 (=)- = (=) «0«0- 4.651 KNm 
17.0930 /2EI 4.8837 /(A4EI 
Me" = +10 + 5082 (=)- = (= ) +0 +0 = 13.663 KNm 








RO 4.8837 (3EI\ , 2.4419 (2EI B 
Me? =—10+0— == (=) E (=) - 0-2 13.663 KNm 


You can compare these answers with answers from force method. 
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